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It is shown that three or more dark/gray sohtons can form bound states in nonlocal media. More 
over dark/gray solitons can form bound states in several balance distances. Numerical simulations 
indicate that some of such bound states are unstable and will decay into a group of fundamental 
solitons, while others may be stable. There exist degenerate bound states with the same velocity, 
Hamiltonian, particle numbers and momentum but decaying in different ways and having different 
lifetimes. 
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I. INTRODUCTION 

As were pointed out previously by N. I. Nicolov, et. 
al.[l[ and Y. V. Kartashov, et. al.l2j, two dark/gray 
solitons can form bound states due to their long-range 
attraction in nonlocal nonlinear media. In this paper 
we will show that three or more dark/gray solitons can 
form bound states also. As were indicated previously [1], 
nonlocal dark/gray solitons have exponentially decaying 
oscillatory tails which in turn give rise to widely dis- 
tributed exponentially decaying oscillatory light-induced 
perturbed refractive index. As a result, in this paper, we 
will find that dark/gray solitons can form bound states 
in several balance distances. Numerical simulation shows 
that fundamental dark/gray solitons are stable. Some 
of such bound states are unstable and will decay into a 
group of fundamental solitons, while others may be sta- 
ble. Bound states of no central symmetry or asymmetry 
will have degenerate states with the same non-vanishing 
velocity, Hamiltonian, particle numbers and momentum. 
Some of such degenerate states are unstable and will de- 
cay in different ways and have different lifetimes. 



II. BOUND STATES OF MULTIPLE 
DARK/GRAY SOLITONS 

The propagation of a paraxial optical beam in a 
medium with a self-defocusing spatially nonlocal nonlin- 
earity can be described by the dimensionless nonlocal 
nonlinear Schrodinger equation(NNLSE) J^-^] 



+ ~ " / - x'\)\u{x', z)|W - (1) 



where u{x, z) is the complex amplitude envelop of the 
light beam, I{x,z) = \u{x,z)\'^ is the light intensity, 
X and z are transverse and longitude coordinates re- 
spectively, R{x) is the real nonlocal response function 
and satisfies the normalization condition / R{x)dx = 1, 
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n{x, z) = — f i?(|a; — x'|)|it(a:', z)\'^dx' is the light-induced 
perturbed refractive index. Note that not stated other- 
wise all integrals in this paper will extend over the whole 
X axis. It is easy to prove if u{x, z) satisfies Eq. ([T]) with 
a nonlocal response function R{x), the scale-transformed 
function u{x,z) — satisfies Eq. ([T]) with an- 

other nonlocal response function R{x) — ^R{^), where 
a is an arbitrary positive real number. So it is adequate 
to consider dark/gray soliton having asymptotic behavior 



\u{x,z)\'^ 



>±oo 



1. Not stated otherwise all dark/gray 



solitons in this paper have this asymptotic behavior. 

By introducing a transformation u{x, z) = tpix, z)e~ 
equation ([T]) turns into 



dt/j 1 d-^ip 



i?(|x-a;'|)|i^(a;')l'dx'- 1 



0(2) 



which has three integrals of motionj: 
number of particles 



namely, the 



(3) 



the momentum 



P = - 
2 



dip* 



dx 

and the Hamiltonian 



■ r— 

dx 



(4) 



H = 



1 



dip 



dx 



dx 



{l-W)Ax 



R{\x-x'\) \ip{x)\^\il:{x')\^ ~1 dxdx'{b) 



In this paper we numerically solve Eq. ^ to find 
dark/gray soliton solutions 



iP{x,z) = (t){x-vz)e'^^''-''''\ 



(6) 



where v is the velocity of the dark/gray soliton relative to 
the cw background, and real functions cp and 9 asymptot- 
ically approach cp{x) °°> 1, 9{x) °°> t9o, where 
6*0 is a real constant. By inserting Eq. ([6|) into ([2]), we 
get 



e' = v[l-^^ (7a) 
^-(pjR{x~ x')(p^ {x')dx' ^ (7b) 
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FIG. 1: (la), (2a), (3a), (4a) are the intensity profiles of 
■i/'OiV'i 1^1,1 iV'i, 1,1 respectively, (lb), (2b), (3b), (4b) are their 
counterpart evolutions. Here w = 5, i; = 0. 
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FIG. 2: Functional dependence of momentum P, Hamiltonian 
H and the number of particles A'' on the velocity v of ipo (solid 
line) and (dashing line) when w = 5. 
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FIG. 3: It needs a longer and longer propagation distance for 
the decaying of bound states ipi as the characteristic nonlocal 
response length decreases from (a) it; = 4 to (b) ui = 2, to (c) 
w = 1.5, and to (d) w = 1.2. 



As an example, we consider this following nonlocal case 
in which the light-induced perturbed refractive index is 
governed by 



which results in n{x,z) = — J R{\x — x'\)\u{x' , z)\'^dx' , 
where i?(|x|) = ^^^V (~''^) ^^'^ ^ characteris- 
tic nonlocal response length of the media. Numerically 
solving Eqs. (O, we can find the dark and gray soli- 
ton solutions. The intensity profiles of dark solitons 
■00 and ^/'i, "01.1: V'l, 1.1 a-nd their counterpart evolutions 
are shown in Fig. ([l} when w = 5,v = 0. (Simi- 
larly, bound states V'li V"!,!; ^i,i,i can form on the in- 
terface of two meida[§]). No observable changes in the 
intensity profile of the fundamental soliton i/jq can be 
found during its propagation . Numerical simulation (not 
shown here) indicate that an initially broadened beam 
Tp{x,0) — tPo{y^,0) will converge into 00 quickly dur- 
ing its propagation. So the fundamental soliton tJjq is 
stable. However, from Fig. ([ij, with no initial pertur- 
bation bound states ■01 1 V'l.ij "01,1,1 are all unstable and 
will decay into a group of fundamental solitons. Numer- 
ical simulations (not shown here) also indicate that ini- 
tially broadened bound states, instead of converging into 
bound solitary states, will ultimately decay into a group 
of fundamental solitons. 

As has been shown previously Q, when w > 1/2 the 



maximal velocity of gray soliton is f max — y ^Xw'^ ■ 
when w = 5, we have Wmax = 0.436. As shown in 
Fig. (I4])(lb), when w = 5 the bound state with an 
initial velocity v = 0.4 is unstable. Numerical simula- 
tions (not shown here) also indicate that "01 with other 
velocities v — 0.1, 0.2, 0.3 are all unstable. So we strongly 
guest bound state 0i with any velocity are unstable when 
w — h. Similarly, bound state 0'i,i, -01,1,1 with any ve- 
locity are all unstable when = 5. When w = 5, the 
functional dependence of momentum P, Hamiltonian H 
and particle numbers iV on the velocity v of 0o and -01 
are shown in Fig. ([21), from which, both for 0o and '0i, we 
get dP/dw < for all velocity. -00 is stable but -01 is not. 
So the criterion of dark soliton instability dP/du > O^J^I 
may be a sufficient but not necessary condition. Namely 
we cannot tell a dark/gray soliton whether stable or not 
if dP/dt; < 0. 

As is shown in Fig. ([3]), it needs a longer and longer 
propagation distance for bound states V^i with w = 
to decay into fundamental solitons as the characteristic 
nonlocal response length w decreases from w = 4 to = 
2, and to w = 1.5, and to w = 1.2. So it is possible 
that bound states 0i could be stable for a small enough 
value of vj. But limited by the nature of the numerical 
simulations method used by this paper we cannot present 
an exact proof for the stability of such bound state. So it 
is still an open question of the stability of bound states 
of dark/gray solitons for small value of w. 
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FIG. 4: (la), (2a), (3a), (4a), (5a), (6a), (7a), (8a) are the in- 
tensity profiles of i/;i,'!/;2,'!/'3,V'4,V'5)V'6,''/'7,V'8 respectively. 
(2b), (3b), (4b), (5b), (6b), (7b), (8b) are tfieir counterpart evo- 
lutions figured in frames moving with velocity v — 0.4 and 
(lb) figured in a frame moving with v = 0.38. Here w — 5 
and velocities of bound states are all v = 0.4. 



III. BALANCE DISTANCES BETWEEN 
DARK/GRAY SOLTIONS 

As was indicated previously Q, due to the nonlocal 
nonlinear response of the media the dark/gray solitons 
have exponentially decaying oscillatory tails when |u| < 

iorw>\. Let cj,{x) 
\x\ — >■ oo, we haveQ 



1 — x{^)^ then when 



xix) « ci exp(-Ai|a;|) cos(A2|a;| + C2) 
where ci and C2 are two constants, and 



Ai 



A2 = 



1 — 'iw'^v'^ + 4i/;Vl — v'^ 



-1 + 4w2w2 ^ - 
4^2 



(9) 

(10) 
(11) 



and 27r/A2 is the oscillatory spatial period. As shown 
in Fig. (I5|)(a), when w ~ 5,v — 0.4, the fundamen- 
tal soliton xo ha-s a serial of maximums and mini- 
mums interdigitally located at Xq = 0,xi = 4.57, a;2 = 

10.29,2:3 = 15.61,2:4 = 21.10,2:5 = 26.52,2:6 = 

31.97,2:7 = 37.41,2:8 = 42.85. Such exponentially de- 
caying oscillatory tails in turn give rise to an exponen- 
tially decaying oscillatory light-induced perturbed re- 
fractive index, and dark/gray solitons can form bound 
states in not only one but several balance distances. 
For example, as shown in Fig. ([4]), when w = 5,v — 
0.4, the balance distances (the distance between two 
deepest dips) of bound states ■01jV'2 7''/'3iV'4jV'5:'06j'07!V'8 
are di = 4.88,10.48,15.97,21.33,26.82,32.24,37.69,43.13, 
where index i runs from 1 to 8. Obviously we have 
di ~ Xi. Let Adi = d^+i — di, we have Adi — 
5.6,5.49,5.36,5.49,5.42,5.45,5.44. On the other hand, 
from Eq. (fTTj) . the half of the oscillatory spatial period 
7r/A2 = 5.44 very closes to Adi. Other cases with differ- 
ent w and V also show the same relation between 7r/A2 
and Adi- To study the interaction of two fundamental 
solitons Vo separated by a distance d, we introduce a 
coupled state ipdix) = (fidix) eTq)[i0d{x)], where 



4'd{x) 



1-Xo (^-i)-Xo i^+l), 
9'=v(l-- 



(12a) 
(12b) 



and the index d of i/jd denotes the distance. Obvi- 
ously, when d — > 00, there will be no overlap between 
Xo(a; — |) and Xo{x + |), and ipd will decouple into two 
well separated fundamental solitons ipQ. We note, when 
d = 2:1, 2:3, X5, • • • , the maximums of Xoi^ — f ) will over- 
lap with the minimums of xo(a^ + f )j and as Fig. (lS])(b) 
shows, the Hamiltonian of 4'd takes the minimums, while 
when d = X2, X4, Xq, ■ ■ ■ , the maximums of Xoi^ ~ ^) will 
overlap with the maximums of Xo(2:+ f ), and the Hamil- 
tonian of tpd takes the maximums. On the other hand, 
as indicated by Fig. (15])(c),(d), the difference between ip^ 
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FIG. 5: (a)fundamental soliton xo{x) (red line) and its max- 
imums and minimums (black dots), (b) the Hamiltonian of 
ipd (red line) and those of tl)3,tp4,i^5,^pe,i^7,ips (black dots). 
(c),(d) the intensity profiles of '!/j2 , 'i/'s (red line) and '!/'d2 > V'da 
(blue line). Here w = 5,v — 0.4. 



and ipda is small, and the difference between ipi and ipdi 
will decrease as di increases. So, as Fig. (©(b) shows, 
though there is a difference between the Hamiltonian of 
ipdi ^n.d those of ipi, the Hamiltonian of ^pi assumes a 
similar behavior of that of ipd, which may qualitatively 
explain why di w Xi and Adi sa Numerical simu- 

lations shown in Fig. Q indicate that "02 , "04 : V'6 , V's are 
all instable, while 03, 05, "07 may be all stable, and "01 is 
weakly instable when w = 5,v — 0.4. 

Numerical simulations with initially broadened bound 
states with vanishing velocity f = 0, seeing Fig. ©, also 
indicate that slightly initial departure from the bound 
states 02 and V'4 can result in a great different evo- 
lutions in the long run. ip2 and '04 are both instable 
when w = 5,v = also. On the other hand the ini- 
tially broadened bound state of ■03 , instead of converging 
into ^03 or decaying into two fundamental solitons, will 
fall in a seemingly eternal (over a distance larger than 
z — 100000) vibrations around 0'3. So bound state ^03 
may be oscillatory-stable. But it needs a more rigorous 
proof beyond the numerical simulation method of this 
paper to judge the stability of 0^3, -05, ipi. 

Since dark/gray solitons can form bound states at sev- 
eral balance distances, we can construct bound states 
like "01.2,1 or ■0^33, and so on. Interestingly, as shown 
in Fig. ([7|), bound states with no central symmetry or 
asymmetry will have two degenerate states, like '0i_2 and 
"02,1, moving with the same non- vanishing velocity and 
having the same momentum, Hamiltonian and the par- 
ticle numbers but decaying in different ways and having 
different life-times. However, there may still exist stable 
degenerate states, like ■03^5 and "05,3 • 



FIG. 6: (la), (2a), (3a) are intensity profiles of %p2,tp3,^4', 
Numerical simulations of initially broadened bound 
states {lh)Mo^),{ic)MTwy' (2b)V'3( ofg ),(2c)^3(t^); 
(3b)'04(Q^),(3c)04(Y^) respectively. Here w = 5,v = 
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FIG. 7: (a) The intensity profile of "02,1; Two degenerate 
bound states (b) ^2,1 and (c) i/)i,2 decay in different ways. 
Here w — 5,v = 0.25. 



IV. CONCLUSION 

In nonlocal media, multiple dark/gray solitons can 
form bound states in several balance distances. Some 
of such bound states are unstable and will decay into a 
group of fundamental solitons, while others may be sta- 
ble. There exist degenerate bound states with the same 
velocity, Hamiltonian, particle numbers and momentum 
but decaying in different ways and having different life- 
times. 
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